Kondo and anti-Kondo resonances in transport through nanoscale devices 
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We study the current through a quantum wire side coupled to a quantum dot, and compare it with 
the case of an embedded dot. The system is modeled by the Anderson Hamiltonian for a linear 
chain, with one atom either coupled to (side-dot) or substituted by (embedded dot) a magnetic 
impurity. For realistic (small) hopping of the dot to the rest of the system, an exact relationship 
between both conductivities holds. We calculate the temperature dependence for moderate values 
of the Coulomb repulsion U using an interpolative perturbative scheme. For sufficiently large U and 
temperature greater than the Kondo temperature, the conductance as a function of gate voltage 
displays two extrema. 
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I. INTRODUCTION 

Transport through quantum .dots have been studied in 
a series of recent experiments,!^!!] in which the peculiar 
features of the spectral density in the Kondo effect are 
manifested. An important fact in these experiments is 
that only one Kondo impurity is present. 

While these experiments correspond to embedded 
quantum dots, very recent studies provide predictions for 
other situations, in particular a quantum dot side cou- 
pled to a quantum wire.QQ In Refs. 4 and 5, decoupling 
approximations were used for infinite Coulomb repulsion 
U, which have limitations at finite temperatures .13 F 
example, in the slave-boson mean-field approximation 
there is an unphysical phase transition at which the im- 
purity decouples from the rest of the system, either when 
the temperature T = Tk (Kondo temperature) or for 
£F _ £d ^ A, where Ef is the Fermi energy, Sd is the im- 
purity level and A = E/2 is the half- width half-maximum 
of the resonant levelJa The method, used by Torio et al. 
is expected to be very accuratejS'EI but it is limited to 
T = 0. 

Here we report calculations of the conductance as a 
function of temperature and gate voltage for embedded 
and quantum dots for U < 8 A. In some.-jexperimental 
situations A ~ 0.15 meV, U ~ 0.6 meV.EKl Thus, this 
work complements previous ones for infinite U or T = 0, 
and presents a more systematic study of the tempera- 
ture dependence for moderate. U. We use second-order 
perturbation theory in c7ll9il3 modified to ensure that 
the approximation is exact in different limits, including 
U +oo.Ert£l 

The paper is organized as follows. In section II we 
write the relevant expressions which relate the conduc- 
tance with the spectral density at the dot. The approxi- 
mation for the latter is briefly reviewed and discussed in 
section III. Section IV contains the main results. Section 
V contains a brief summary and discussion. 



II. CONDUCTANCE FOR EMBEDDED AND 
SIDE DOTS 

For both cases, the Hamiltonian can be written as an 
Anderson model 



H — H + H con , 



(1) 



where Hq describes a tight-binding chain without site 
plus a disconnected quantum dot, and H con contains the 
hopping terms which connect the quantum dot with the 
chain, 



H = -t ( C iad-la + H - c ) + £ dYl d Ua 



-Ud^d-[d\di. 



(2) 



For the embedded dot, 

Kon = ~ J2 [4 + t[c la ) + H.C.] , (3) 

while for the side dot, 

Kon = ~ E [ C 0^ (*-l C "l* + ^Cla + t'd a ) + H.c] . (4) 
a 

P^he conductance for an embedded dot can be written 

df 



du> 



r CT (o;)p (2<7 (a;), 



(5) 



where f{u>) is the Fermi function, pdo{u) is the spectral 
density of the d a states, and 

r i \ r_i CT (w)ri CT (w) . , 2 , . 

r<rM = — — — —, with T i(J {uj) = 27r(iJ p icr {iv) 



r_i CT (w) + r la (uj) ' 



(0) 



where pi a {uo) is the spectral density of the states de- 
scribed by the operators Ci a in Hq. From the equations 



1 



of motion for a semi-infinite chain, one easily obtains 

1/2 

P-ier = pier = [l — (ui/2t) 2 ~\ /irt. The resonant level 
width of the d states is A CT (w) = 7r^ (t^) 2 Pioi^)- In the 

i 

following we assume that t[ is small, implying that the 
Kondo temperature is much smaller than the band width. 
The sum over spin indices reduces to a factor 2, since for 
zero magnetic field there is no spin dependence. Then, 



G, 



2e 2 47rA e (t[t'_x) 
~[W) 2 + (i'-i) 



£j p d », (7) 



with 



(*i) 2 + y_i) 2 
t 



1 - 



(8) 



where ep is the Fermi energy. A e is the half- width at 
half- maximum of pda(u) for U = 0. Both A e and t/ are 
the parameters of the perturbation calculation, which al- 
lows to obtain pd a (to) as described in the next section. 

For the side dot, Eqs. (||,[?]) apply with the substitution 
Pda — * Poa, and t\ — > t. Assuming t-\ = t\ = t (periodic 
chain) for simplicity and t' <C t as before, we have 



G„ = 



2e 2 



hp(eF) 



df \ 



(9) 



where p(u>) = 1/ 



7T Ut 2 



,1/2' 



is the density of states 



of any site of the chain for given spin and t' = 0. The cor- 
responding density of states at site for t' ^ (poo-(w)) 
can be expressed in terms of pdo-(u>). Using equations of 
motion for the Green's functions, one can write: 

G 0a (f, w) = G O a(0, u) + (t') 2 G d a(oj)Gl a (u), (10) 

where Go a (t' , w) = ((co a ; Cq„))u is the Green function for 
site and hopping t with the dot. Similarly Gda{t' , to) = 
{(d a ; dl)) u is the Green function for the electrons in the 
dot. From Eq. ( |l0| ) and Goct(0,w) = — i7rp(u;),one has 

1 2 
Po<tM = ImG 0(r (t,w) = - [7rfp(u>)] pdo-M- 

7T 

(11) 

Replacing Eq. (||) and neglecting again the dependence 
of p(to) on w near e^,we can write 



G 



1tA s J Pda(u) 



where 



A. = 



{t'f 



y^t 2 



(12) 



(13) 



We see that in both cases, G e and G s are determined 
by an integral of the same form, involving the density 
of states at the dot pd a {^>), but with different sign. An 



increase in pda(w) for an embedded (side) dot leads to an 
increase (decrease) in G e (G s ) . 

At T = 0, the conductivities can be related with the 
occupation of the .impurity level nd a — {d\,d a ) , using 
Friedel's sum ruleO For e e ff and A independent of fre- 
quency, this rule states that 



Pda{£F) = 



sin 2 (7rn fiCT ) 
ttA 



(14) 



and then, replacing in (^) and (|T^) one obtains the simple 
result, 

_ 2e 2 4(t^Q 2 . 2 

G e = —r-- 72 sin (Tm dcr ), 



n 2e2 K 

G s = —— COS {TTTldo 

h 



(15) 



Thus, the ideal conductance for an embedded dot is ob- 
tained for a symmetric dot t'_ 1 = t[ and when V g is tuned 
to the symmetric Anderson model (e e ff = Ep). For the 
side dO|t~-the dependence with V g is the opposite, as noted 
earlier .Ou 



III. THE SPECTRAL DENSITY 

The dot Green's function can be written in the form 

1 



G dl r(uj) 



to - Eeff +iA- £<r(k>,T) 

For the embedded dot, A is given by Eq. (@) and 



(16) 



e% ff =s d - eV g + Un^-^ [(t[) 2 + (t'^) 2 ] , (17) 

where the subscript a in £^ff nas t» een dropped for sim- 
plicity and V g is a gate voltage which controls the dot 
level. The dependence on frequency of the terms pro- 
portional to (t-) 2 has been neglected assuming t' <C t as 



before. For the side dot, A was written in Eq. (pL3|) and 



'-eff 



£d ~ eV„ + Un, 



(1 n . 



(18) 



In traditional second-order perturbation theory in U, 
the self-energy H a (uj,T) is calculated from a Feyn- 
mann diagr.arn_which involves two sums in Matsubara 
frequencies:Ejil3 

4 2 >(iw,,T) = -(UT) 2 Y,G Q dA^i-^m)GU^n) 

xG%(w n + iv m ), (19) 
where the unperturbed Green's function is: 



G da (io) = 



UJ — £ 



eff 



■iA 



(20) 
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The choice of * is arbitrary and is equivalent to the 
division of the flamiltonian between unperturbed part 
and perturbation. Taking = £ e // [ js^quivalent to a 



sum of an infinite series of diagrams .EJE3 If the expec- 
tation value nda were calculated using G da , e e jf would 
represent the Hartree-Fock effective d level. However, we 
find that in general calculating-p,^ with Gd a {^>) agrees 
better with Friedel's sum ruleJlH Then 



nda 



f(u>)pda(u)du, with pd a (oj) = ImG^fu). 

7T 



rida should be determined self-consistently from Eqs. 
0.0 . and (1^) or (|J), for each value of V g . 
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perature dependence of the Kondo peak for T < 2Tk 
agrees within 20% with WRG results. The temperature 
dependence of the rest of pdo is very small. 

As e e f f increases from the symmetric situation e e f / = 
Ep, the Kondo peak is slightly shifted to higher ener- 
gies, but Pda(£F) remains high. However, for U > 5A, 
the perturbative result shifts to the opposite direction, 
leading to a violation of Friedel's sum rule. This short- 
coming can bexuped using an interpolative solution for 
the self-energy£ 3 



(21) E ff ( Wj T) 



(2) 

rida(l - ndwjT^ 



)-[{l-n dw )U + e d r 



-2v( 2 ) 



(22) 



where 



= — J f(u>)Im.G d \ 7 (u))du>/TT is the impurity 
occupation corresponding to the unperturbed impurity 
spectral density, and e\ — e e // — Uridw is the effective 
impurity level for U — (see Eqs. (0) and (|lq)), which 
depends on the applied gate voltage V g . This self energy 
leads to the exact Gda(u) not only for U — 0, but also 
for a decoupled dot (A = 0), and reproduces the leading 
term for lo — > oo. The effective unperturbed impurity 
level e° e ff can be calculated selfconsistently at arbitrary 
temperatures imposing the condition n d(T = n^EjEJ. 
Unless otherwise stated, the results presented here are 
obtained following this approach. For T = 0, the results 
can be improved further determining Egtt by imposing 
Friedel's sum rule, as done by Kajueter and Kotliar to 
solve the impurity problem in the dynamical mean-field 
approach.Ej We call T,kk the self energy obtained in this 
way. 
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FIG. 1. Spectral density at the dot for given spin, U = 4-7T, 
Ed = — 2-7T, ef = 0, and several temperatures. A = 1 is chosen 
as the unit of energy. 

A particular situation is the case e e ff = Ep, at which 
Pdcr(w) is symmetric around ep and then nd a = 1/2. 
At this point, the-perturbative result satisfies exactly 
Friedel's sum ruleO Out of this point, for large U, the 
perturbative result should be improved as discussed later 
in this section. For the symmetric model {s e ff — £f), 
the spectral density coincides with that obtained using 
Quantum Monte Carlo with Maximum Entropy method 
within statistical errors if U/A < l.25n ~ 3.93.N For 
U I A ~ 7.6, the difference between both results is ~ 10%. 
In Fig. 1 we show the spectral density in the symmet- 
ric case for U/A = 4ir. Although this ratio is beyond 
the validity of the perturbation theory, the result is in 
qualitative agreement with the corresponding result oh= 
tained using Wilson's Renormalization Group (WRG),Ej 
The central (or Kondo) peak is broader in perturbation 
theory. However, if the Kondo temperature is defined 
as half width at half maximum (Tr- ~ 0.10A), the tem- 
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FIG. 2. 

In Fig. 2 we show the spectral density for U = 8A 
in an asymmetric case with E d — Ep = —3 A. The spec- 
tral density shows the characteristic charge fluctuation 
peaks near Ed and Ed + U, and the Kondo peak near 
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£p. The slave-boson mean- field theoryi cannot repro- 
duce the charge fluctuation peaks and brings an incom- 
plete description of the spectral density. For the sym- 
metric case Ed — £f = —4 A, the Kondo temperature de- 
fined as half-wjxilh at half-maximum of the Kondo peak 
is Tk — 0.23A.E3 In a scale of Tk, for increasing temper- 
atures, the Kondo peak is rapidly suppressed, while the 
charge fluctuation peaks remain and absorb the spectral 
weight of the Kondo peak. For T ~ Tk, the density of 
states at the Fermi level Pda(£F) is reduced to approxi- 
mately half its value at T = 0. 




FIG. 3. Dashed line: dot occupation for a given spin as 
a function of gate voltage for (7 = 8, and T = 0. Dotted 
line: corresponding result obtained fromr-the phase shift at 
the Fermi level using Friedel's sum rule.E-j Full line: result 
using Ekk (imposing Friedel's sum rule). The zero of V g is 
set at the point where e e // = ef, rida = 1/2. 

In Fig. 3 we compare the resulting-p^er at T — 0, with 
that obtained from the phase shiftO for the largest U 
(worst case) considered in this work (U/A = 8). We see 
that the deviations are small and concentrated near but 
out of the symmetric case. Note that the result using 
Tikk (in which Friedel's sum rule is imposed) is very 
similar to that obtained imposing n da = rider- In turn, 
this result agrees r^fry well with exact calculations using 
the Bethe ansatzJiij 



IV. RESULTS FOR THE CONDUCTANCE 



From Eqs. (§), (§), ©, @, ©, (O), @, and 
(pi|), one can see that for given U and t' (or i-) <C t, 
the conductance for the embedded and side dots depend 
on two parameters: A and Ed — eV g . Furthermore, in- 
dependently of the approximation used for the spectral 
density, both conductances are related shifting V g (or Ed) 
and rescaling A : 



G (A V) = mt '- 1 ' 
y ' B> [(t[) 2 + (<'_i) 2 ] 2 



2e^ 
h 



G s (A,V g ) 



, (23) 



with 



v; = K, + [(t' 1 ) 2 + (t'_ 1 ) s 



■g -y • L v-i, , v--i, j 2t2e - 

Thus, it is sufficient to calculate one of the conductances. 
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FIG. 4. Conductance for the side dot in units of 2e 2 /h as a 
function of gate voltage for several values of U and T indicated 
inside each figure. The bottom curve with thiner full line was 
obtained using Y>kk- 

In Fig. 4 we represent G s for several values of U 
and T. The unit of energy is set as A = 1, the con- 
ductance is expressed in units of 2e 2 /h and the zero of 
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the gate voltage is set at the symmetric Anderson model 
(eVgo = e d + U/2 = for the side dot). For T = we haj^ 
also included the result imposing Friedel's sum rule,t3 
(using Yjkk) which gives a slightly smaller value of G s 
and allows to estimate the error at finite T. This error is 
due to the effect on pdai^F) in the real part of £(ef) for 
the asymmetric Anderson model. The imaginary part is 
always zero as it should be for a Fermi liquid. For U = 8, 
the difference in G s at T = between both interpolative 
methods attains its maximum value ~ 0.07 x 2e 2 /h near 
eV g — eVgo = ±2.5. For U = 5, the maximum deviations 
are of the order of 2%, and for smaller U they are negligi- 
ble. Note that the difference between both approaches is 
much smaller in integral properties, like n^r (see Fig. 3). 
As a consequence G s obtained using Eq. (|lq ) (instead of 
Eq. (|l2|)) is practically the same in both approaches. Us- 
ing T^kk, of course Eq. ([l^) and Eq. ([ll]) give identical 
results. 

The conductance for T = and U > 5 are qualita- 
tively similar to the corresponding result of Ref. 6. As 
the temperature is increased, G s near V g — rapidly 
increases. This is expected from the temperature de- 
pendence of the Kondo peak, presented in the previous 
section. For U > 5, at temperatures of the order of the 
Kondo temperature [T K ~ 0.45 for U = 5 and T K 0.23 
for U — 8 defined as half width at half maximum of the 
Kondo spectral peak), or slightly smaller, the conduc- 
tance shows two minima. For T > Tk, the resulting 
structure is qualitatively what one expects for t' — > 0, in 
which only the Coulomb blockade peaks at £4 and Ed + U 
are important at finite temperatures. However, the peaks 
are displaced towards zero, particularly for small U and 
T, For U = 5 and T > Tk, the interpolative scheme 
and the ordinary second order perturbative result (using 
£<j = Sct and with sclfconsistency in ride only) are prac- 
tically identical. At smaller temperatures, the ordinary 
treatment exaggerates the structure with two minima. 
Instead, the comparison with the results using ^kk at 
T = and [7 = 8, suggest that the interpolative treat- 
ment with rida — n° da inhibits somewhat this structure 
for large U. 




FIG. 5. Temperature dependence of the conductance for 
the embedded dot and the density of states at the Fermi level 
for the symmetric Anderson model and U = 5A. 

As U is decreased, the structure with two minima 
weakens. For U = 3 A the Anderson model is near 
the intermediate valence regime, since the Kondo regime 
requires U A. While the spectral density displays 
only one peak, without marked shoulders for [7 = 3, 
plateaus reminiscent to the transition to the Coulomb 
blockade regime are still present in the conductance for 
0.2 < T < 0.5. In Fig. 5 we show the temperature depen- 
dence of the density of states at £f and the conductance 
of the embedded dot for V g — and [7 = 5. The shape 
of both quantities in a logarithmic temperature scale is 
similar to the experimentally observed conductance.0 In 
our case the linear behavior with logT extends over less 
than an order of magnitude: from T ~ 0.13 to ~ 0.5 for 
Pdai^F) and from ~ 0.1 to 0.3 for G e . This corresponds 
to values lower or of the order of Tk- 
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FIG. 6. Same as Fig. 6 for U = 3A. 
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In Fig. 6, we show the same temperature dependence 
for U — 3. In this case, the half width at half maximum of 
the spectral density in the symmetric case is 0.73. In com- 
parison with Fig. 5, the behavior of Pda(£F) and G e vs. T 
is qualitatively different. pd a ( £ F) ~ log? 1 for 0.2 < T < 
0.6. In contrast, except for a very small curvature up- 
wards, G e is linear in logT for 0.16 < T < 2. A quadratic 
fit of all points with temperatures an integer times 0.01 
in this interval gives G e h/(2e 2 ) = 0.4295 - 0.2112 log(T 
/A) + 0.0050[log(T /A)] 2 , with a mean square deviation 
6 x 10~ 5 . 



V. SUMMARY AND DISCUSSION 

We have calculated the conductance for a quantum dot 
embedded in a quantum wire or side coupled to it as a 
function of gate voltage and temperature.. We have con- 
sidered moderate values of U which have not been studied 
before. For these values of U, the perturbative method 
used is quite accurate. For small coupling of the dot to 
the wire, the conductance is determined by an integral 
of the density of states at the dot times the derivative of 
the Fermi function. As the temperature is increased, for 
U > 4A we find that the conductance as a function of 
the gate voltage changes from a typical shape with one 
extremum dominated by the Kondo peak in the spectral 
density, to another with two extrema, corresponding to 
the Coulomb blockade regime. To our knowledge, this 
crossover has not been reported previously. 

For U = 3A, we find that the conductance can be 
described as linear in log T as a very good approximation 
over more than an order of magnitude in T. This behavior 
is rather unexpected and difficult to explain in simple 
terms, since different energy scales in the problem are of 
the same magnitude. 
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